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Abstract 

It has recently been observed that the weakly coupled plane wave matrix 
model has a density of states which grows exponentially at high energy. This 
implies that the model has a phase transition. The transition appears to be 
of first order. However, its exact nature is sensitive to interactions. In this 
paper, we analyze the effect of interactions by computing the relevant parts of 
the effective potential for the Polyakov loop operator in the finite temperature 
plane-wave matrix model to three loop order. We show that the phase transi- 
tion is indeed of first order. We also compute the correction to the Hagedorn 
temperature to order two loops. 
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1 Introduction 



It is conjectured that the full dynamics of M-theory is encoded in a certain 
maximally supersymmetric matrix quantum mechanics called the BFSS matrix 
model 1 . At finite N that model is thought to describe a discrete light-cone 
quantization of M-theory in the infinite momentum frame. The light-cone is 
decompactified by a taking a particular large N limit where the discrete light- 
cone momentum p + = N/R is held constant, with R the radius of the null 
identification x~ ~ x~ + 2ttR. 

Recently, a massive deformation of the BFSS matrix model appropriate to 
M-theory on a plane wave background has been studied^. The main differ- 
ence between it and the BFSS model is their super symmetry, which in the 
latter case is that of a maximally supersymmetric plane wave spacetime rather 
than 11-dimensional Minkowski space. A great advantage of the plane wave 
matrix model is that it has a weak coupling regime where it can be studied 
systematically using perturbation theory |3*]-|11|. It also has a powerful super- 
symmetry algebra which allows the extrapolation of some perturbative results 
to the string coupling regime. 

The action of the plane wave matrix model is 

S = — / drTr DX i DX i + DX~ a DX~ a + i^ Ia D^ Ia + — [X\ X j ] 2 + 
2R J \ 2 

+^-[X a , X 1 } 2 + R 2 [X a , X 1 ] 2 + ^[X\X j ] 2 + ^-[X a , X 1 ] 2 + R 2 [X a ,X' 1 } 2 
-R^ Ia a^[X\^) + ^e a ^ aI ^ u [X\^ J ] - ^e a ^ aI {^) IJ [X\^ aJ ] 



f ) 2 {x ~ a)2 - (I) 2 {xi)2 - 1^ 1 * - ^ 2 4 e ^ x ~ axbx5 ) (i) 



where all of the variables are N x N matrices. The indices on the bosonic 
matrices take values a,b,c = 1,2,3, i,j = 4, . . . , 9. The U(N) symmetry is 
gauged. All variables transform in the adjoint representation of the gauge group 
X % — > UX i W, etc. The time derivatives are covariant, D = d T — i [A, ...] with 
an N x N Hermitian gauge field A. The fermions have 8 complex components 
with I, J = 1, . . . , 4 and a,/3= 1,2. The spin matrix has the property g J (g' ? )^ + 
g^g*)^ = 25^14x4. e a p and are antisymmetric tensors. 

Note that the last line of Q contains terms with the parameter /i which 
comes from the plane wave geometry. That geometry has metric and constant 
four- form, 

ds 2 = -2dx + dx- + dx a dx a + dx l dx l - dx + dx + {x a ) 2 + (|) 2 (x*) 2 
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It is a maximally symmetric solution of 11-dimensional supergravity for any 
value of the parameter ju. If we set /i to zero we recover the Minkowski metric. 

The plane wave matrix model Q can be systematically analyzed in pertur- 
bation theory. When N is large, the expansion parameter is ^2] 



Note that all dimensional quantities are in units of the eleven dimensional 
Planck length which we have set to one. 

Aside from being a model of M-theory on a plane-wave background, there is 
an interesting connection between the matrix model (J2) and four dimensional 
J\f = 4 supersymmetric Yang-Mills theory. At the classical level, the degrees 
freedom and dynamics described by (^Q) are a consistent truncation of the Yang- 
Mills theory which keeps only those modes which are invariant under a certain 
SU(2) subalgebra of the full SU(2, 2|4) superconformal algebra. In the planar 
limit, this truncation seems to also hold at the one- loop level (§]. It is then 
natural to speculate about whether it would inherit some of the properties of 
J\f = 4 supersymmetric Yang-Mills theory, such as integrability of the large 
N limit, at least in the leading orders of perturbation theory. This has been 
studied in detail in recent work It has been shown that the form of 

the dilatation operator for supersymmetric Yang-Mills theory which has been 
conjectured using integrability as an input, when restricted to the appropriate 
sector of the theory, is identical to the Hamiltonian of the plane wave matrix 
model, at least for computing energy levels to three loop order. 

1.1 Thermodynamics 

Some aspects of the thermodynamics of the plane wave matrix model have 
been studied before |13j-jl8j. It was shown in ^3] that, at weak coupling, and 
in the large N limit, there is a deconfinement phase transition at a particular 
temperature. The use of the word "deconfinement" in a theory where there is no 
spatial extent over which particles can be separated must be justified carefully. 
It is possible to distinguish the two phases by the behavior of the free energy, 
F[T]. In a weakly coupled theory, the free energy should be proportional to the 
number of degrees of freedom. In the confining phase the number of degrees 
of freedom at a given energy are of order the number of color singlets, which 
in the lower energy parts of the spectrum does not grow with the rank of the 
gauge group, N, but is of order one. In a deconfined phase, the number of 




(2) 
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degrees of freedom is the number of elements of the matrices, which is of order 
iV 2 . Thus, we would expect 

F 

lim — k = confined 

N^oo N 2 

F 

lim —77 ^ deconfined (3) 

JV^oo N 2 

When interactions are turned on, the large N limit which we are using is the 
't Hooft limit where the coupling constants are also scaled as N is taken 
large so that A defined in (J2J) is held constant. Note that this limit is different 
from the large N limit with N/R held constant which would decompactify the 
null direction in M-theory. 

In a gauge theory such as (|T|). where all variables transform in the adjoint 
representation, there is an order parameter for confinement called the Polyakov 
loop |19[ 120] . It is the trace of the holonomy of the gauge field around the finite 
temperature Euclidean time circle, 

P=±Tr(e i * A ) (4) 

It is related to the difference of free energies of the state where an additional fun- 
damental representation quark is introduced and the state without the quark, 

F q [T] - F [T] = -T ln(P) 

The free energy of a free quark in the confined phase should be large compared 
to that in the deconfined phase. In fact, 

(P) = confined 

(5) 

(P) ^ deconfined 

is interpreted as requiring infinite energy to insert a fundamental representation 
quark into the system when it is in the confining phase, whereas it is finite in 
the deconfined phase. The Polyakov loop operator is widely used to study 
finite temperature de-confining phase transitions in higher dimensional gauge 
theories, 

A phase transition where both criteria for deconfinement @ and (J5J) occur 
has been found in the large N limit of the matrix model ^3]. It has been 
argued to be generic to the large N limit of matrix quantum mechanics when 
the matrices have a gapped spectrum. The transition temperature is of order 
of the mass gap. 

An intuition for why this phase transition occurs can be gained from count- 
ing matrix degrees of freedom As a simple example, consider a set of d 
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matrix harmonic oscillators with Hamiltonian 

d 



H = y^uTr (a\ai 



=1 



and matrix-valued creation and annihilation operators with algebra 

S ij c 



(6) 



(7) 



States are created the creation operators a^ b operating on a vacuum |0 >. To 
get a state with energy E = nuj we must act with n creation operators. 

The analog of gauge symmetry is to require a physical state condition of 
invariance under the unitary transformation 



o 



a 



uau 



where u E U(N). We assume that the vacuum state is invariant. Then, phys- 
ical states are created by operating with invariant combinations of creation 
operators. In the limit N — > oo all such combinations are traces 



Tr (aj 



»i 



Tr ( at at 



n 32 



"2 



Tr [a\a\ 2 al 3 



10 > 



(8) 



where the energy is 

E = uj{ni + 2n 2 + 3n 3 + . . .) 

The number of these traces with a fixed energy, E, does not scale like TV 2 
as N — > oo, instead it approaches a constant as N is taken large. Thus at 
low enough temperatures, the free energy should approach an A^-independent 
constant as is taken large. 

However, the number of independent traces does increase rapidly with the 
energy E. It has been shown ^5] that, in the large limit, the oscillator has 
a Hagedorn density of states at high energy, 

1 



P{E) 

where the Hagedorn temperature is 



,E/T B 



E 



Th = ojJ In d 



(9) 



This mirrors a similar discussion for weakly coupled Yang- Mills theory |31|-|36|. 

At temperatures higher than T#, the thermodynamic canonical ensemble 
does not exist. It could be made to exist by keeping A^ large but finite. That 
would cut off the exponential growth in the asymptotic density of states at 
some large energy, of order N 2 . Then we could consider a temperature that is 
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greater than T#. Both the energy and entropy would be dominated by states 
at and above the cutoff scale and would be of order N 2 . The divergence of the 
free energy ~ iV 2 occurs as we take the limit N — > oo at constant temperature 
(noting that the Hagedorn temperature does not depend on N). 

It was shown in ref. [T^ that, from a different point of view, this behavior 
of the matrix model can also be found as a large N Gross- Witten type of phase 
transition |37l I38j which is familiar in unitary matrix models. In the next 
Section, we will discuss this in more detail. 

2 One loop 

Now, let us consider a perturbative expansion of the plane wave matrix model 
at finite temperature. The partition function has Euclidean path integral rep- 
resentation similar to those of finite temperature quantum field theories. We 
consider the theory in Euclidean space where time is compact and identified as 

r ~ t + P 

with 

P = l/T 

T is the temperature. The partition function is given by the functional integral 

where L is the Euclidean time Lagrangian 
L = — Tr (bX'DX 1 + DX a DX a - ^ Ia Di> Ia 

+ ^ TV ((f) 2 + O 2 {Xlf + I^ Ia + ^4 e ^ X " XbX5 

+ R^ Ia af[X a ,^] - ^e a ^ aI g\j[X\^ J ] + |e^Vw(g lt ) /J [*\ Vw] 

-^[X\xi] 2 - ^-{X a ,X~ b ] 2 -i? 2 [X s ,JT] 2 

(11) 

The bosonic and fermionic variables have periodic and antiperiodic boundary 
conditions, respectively 

A(r + f3)=A(r) , X* {t + (3) = X' (r) , ^(r + (3) = -^(r) 



Since the boundary conditions for fermions and bosons are different, super- 
symmetry is broken explicitly. Of course this is expected at finite temperature 
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where bosons and fermions have different thermal distributions. Supersymme- 
try is restored in the zero temperature limit. We will see the results of this 
explicitly in the following. A parallel discussion of the BFSS matrix model at 
finite temperature can be found in ref. |39j-|4Uj. 



2.1 Gauge fixing 

To begin, we must fix the gauge. It is most convenient to use the gauge freedom 
to make the variable A static and diagonal, 

~rA a \ } = , A a }j = A a 5 a b 
or 

Once this is done, the remaining degrees of freedom of A are the time-independent 
diagonal components, A a . We shall see that they eventually appear in the form 
exp (if3A a ). 

The Faddeev-Popov determinant for the first of these gauge fixings is 1 

(-s H + > (A ° - *>)) = det ' H) det ' Hr + ' (a ° - *>) 

(12) 

where the boundary conditions are periodic with period (5. The prime means 
that the zero mode of time derivative operating on periodic functions is omitted 
from the determinant. Once the gauge field is time-independent, we do the 
further gauge fixing which makes it diagonal. The Faddeev-Popov determinant 
for diagonalizing it is the familiar Vandermonde determinant, 

\{\Aa-A b \ 

This is also just the factor that the time independent zero mode would con- 
tribute to the second of the determinants in 1)12(1 . Including it gives the deter- 
minant 

n det/ (-^) det (-^+^-^)) ( i3 ) 

where there is now no prime on the second factor. These determinants can be 
found explicitly. We will do this shortly. 



^^Using zeta-function regularization, 
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2.2 Classical ground states 

We will perform a semiclassical expansion of the free energy in the sector of 
the theory with the classical ground state X a = = X 1 in the large N, 't 
Hooft limit. This is a double expansion. First, it keeps only planar Feynman 
diagrams, which is the leading order in an expansion in 1/N 2 . Secondly, it 
is perturbative in that it keeps those diagrams which are of low orders in the 
coupling constant. Since N is large, the appropriate coupling constant is A 
defined in eqn. ©. 

The configuration X a = = X % has zero classical energy. In fact, because 
of supersymmetry, the ground state energy of the theory quantized beginning 
with this vacuum is zero to all orders in perturbation theory. This will provide 
us with an important check of our finite temperature computations. Of course, 
temperature breaks supersymmetry and the free energy of the thermodynamic 
state, which is what we compute, is non-zero. However, we will always be 
able to check whether the zero temperature limit of the free energy, which is 
the ground state energy, vanishes. We find that its expansion in the coupling 
constant A indeed does so to order A 2 . 

In this subsection, for completeness, we comment on the fact that there are 
a large number (— > oo as N — > oo) of other ground states with zero energy. To 
see this, observe that the classical potential in (^Q) can be written in the form 



V = -Tr 
2 



(j^X a + ie abc X b X c ) +-(i{X l ,Xi}) 2 + 



+ (i[x\x«]) 2 + ^)\x 

It has isolated global classical minima (V = 0) where 



i\2 



(14) 



Ki = , *ci = ^</ a (15) 

and J a form an N-dimensional representation of the SU(2) algebra, [J a ,J b ] = 
ie abc J c . In addition, the classical solution for gauge field must obey the equation 

[^ci,J c a i]=0 

If J a are an irreducible representation of SU(2), by Schur's Lemma, A c \ must 
be proportional to the unit matrix A& = c ■ I and a symmetry of the theory 
allows us to set c = 0. The gauge symmetry is realized by the Higgs mecha- 
nism and fluctuations of the gauge field are all massive. On the other hand, 
when the representation is reducible, there are gauge fields which commute 
with the condensate. The parts of A c \ which commute with the condensate 
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remain undetermined by the classical equations. The volume of the space of 
all possible such A c \ forms a moduli space of the classical solutions which must 
still be integrated over, even to obtain the leading order in the semi-classical 
approximation to the partition function. 

The implications of the solutions (JT5|) have been discussed in detail in ref. [H]. 
They were interpreted in terms of the spherical membrane and transverse spher- 
ical 5-branes which exist on the 11-dimensional plane wave background. 

In refs. |13[ I17j it was argued that the phase transition that we shall study 
here occurs only when there is a residual gauge symmetry and only in the limit 
where the rank of the residual gauge group goes to infinity. According to ref. |5j , 
this is the limit of the theory which describes 5-branes. Perturbation theory in 
that limit is governed by a 't Hooft coupling similar to © with N replaced by 
the rank of the residual gauge group. Ref. jH] argued that, in the 't Hooft limit, 
the barrier between the degenerate vacua becomes infinitely high. They also 
argued that the limit decouples the 5-brane from other degrees of freedom and 
focuses on its internal dynamics. 

The results of refs. |13l 117] can be interpreted as saying that the phase tran- 
sition of the matrix model occurs only in 5-brane states and not in membrane 
states, and it seems to be associated with internal dynamics of the 5-branes. Of 
course, refs. |131 117j analyzed only the weak coupling limit of the matrix model 
which is far from the limits which are conjectured to describe supergravity of 
an 11-dimensional spacetime continuum where the 5-brane would live. In order 
to apply it directly to any known behavior of 5-branes, it would have to be 
extrapolated to strong 't Hooft coupling. This is of course a difficult problem. 
The perturbative expansion in the present paper perhaps gives an indication 
that the phase transition is of first order but much more would have to be 
done to answer the question of persistence or nature of the phase transition in 
the 5-brane regime. One interesting extension of the present work would be to 
examine the dependence of the nature of the phase transition on the number 
of 5-branes. For k coincident 5-branes, X a contains N/k /c-dimensional repre- 
sentations of SU(2), with k being held fixed in the large N limit, ^-dependence 
of the phase transition temperature was investigated in ref. |13| . 

2.3 Semiclassical expansion 

If we expand about the classical vacuum = = X^, we find the partition 
function in the 1-loop approximation is 




(16) 
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where D a ^ = -4- — i(A a — Ab). The first two terms in the numerator are the 
Faddeev-Popov determinant. The third term comes from fermions whereas the 
denominator is from bosons. Using the formula 

, ( d \ n ■ i 

det h lu = 2 smh — 

V dr J 2 

with periodic boundary conditions and 

. d \ (3uj 
det ( — - — h uj = 2 cosh — 
dr J 2 

with antiperiodic boundary conditions, we can write 2 

N 

- a 

2~7T -LI [1 _ e - / 3/V3+i/3(A a -A i) )]3[i _ £-,3^/6+1/3^ -A b )] 



|-r d(/3A,) TT [1 " e^ A "-^)][l + e -/^/4+i/3(A a -A b )]8 



Note that, because of supersymmetry, the zero temperature (/? — > oo) limit 
of the partition function is one. It also has a symmetry under replacing e - ^ 
by 

We must now do the remaining integral when N — > oo. There are N inte- 
gration variables A a and the action, which is the logarithm of the integrand is 
generically of order N 2 which is large in the large N limit. For this reason, the 
integral can be done by saddle point integration. This amounts to finding the 
configuration of the variables A a which minimize the effective action: 

S eS =J2{~ ln tl " e*^-^] - 81n[l + e -/W4+i/3(A a -A fe )j + 

a^b 

+ 31n[l - e-^/3+^^a-A,)] + 61n[1 _ e -^/ 6 +^(A a -A,)A (18) 



To study the minima, it is illuminating to Taylor expand the logarithms in the 
phases (this requires some assumptions of convergence for the first log) 



^ 1 - 8(-) n+1 r 3n - 3r 4n - 6r 2n . . 

S eS = } <?-n<?n (19) 



71=1 



Here, 



and 



exp (-/3/i/12) 

e in0A a ^ 



1 N 



N ^ 

a=l 



2 Because the matrix model action {Q is invariant under replacing A by A plus a constant times 
the unit matrix, we see that the integrand in eqn. (|17|l is indeed invariant under translating all values 
of A a by the same constant. 
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Recalling (jlj), we note that (j) n are multiply wound Polyakov loop operators 
evaluated in the static, diagonal gauge. The zeroth moment is normalized 

0o = 1 (21) 
The other elements are constrained by sum rules. The density defined by 

-. oo 

p(x) = »?I>(x-/340 ^° 

a=l 

n 

is a non-negative function. For example, if only 4>q and 4>±\ are nonzero, (|22|) 
implies that \(j)i\ < 1/2. 

In this one-loop approximation, the action is quadratic in the Polyakov 
loops. When all coefficients of the quadratic terms are positive, the action 
is minimized by <j) n = for n ^ 0. This is the confining phase. When a 
coefficient becomes negative, the effective action is minimized with one of the 
loops nonzero. The result is a condensation of the loops. 

As we raise the temperature from zero (and lower (5 from infinity), the first 
mode to condense is n = 1. This occurs when 

r c = 1/3 - T c = — f- « .0758533^ 
12 m 3 

and 4>i ^ when T > Tq- There are some speculations about the nature of 
the distribution of the angles (3A a in the deconfined phase in ref. |17| . 

Note that this condensation breaks a U(l) symmetry. This is associated 
with the center of the gauge group U(l) £ U{N). It arises from the fact that 
all variables are in the adjoint representation. In the Euclidean path integral, 
gauge transformations X{t) — > U{t)X(t)U* (r) must preserve the periodicity 
of the dynamical variables. They therefore must be periodic up to an element 
of the center, U((3) = e lS U(0). The Polyakov loop, on the other hand, being 
the holonomy on the time circle, does transform as P — > e ld P. 

Even once the static, diagonal gauge is fixed, there is a vestige of this 
symmetry where j3A a — > (3A a + 9 or <p n —>■ e in6 (j) n . This symmetry restricts the 
form of the effective action for Polyakov loops, so that the term with (j)^ . . . (pk n 
must have ^ fcj = 0. It is a good symmetry of the confined phase and it is 
spontaneously broken in the deconfined phase. The Polyakov loop operator is 
an order parameter for this symmetry breaking. 
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3 Higher loop order 

In the previous Section, we have computed the effective action as a function 
of the variables cj) n in the one-loop approximation. The result is the quadratic 
potential in eqn. (|19|). To find the stable phase, it is necessary to find those 
values of (f> n which minimize (|19|) . subject to the constraint (|22j) . 

The phase transition occurs when the curvature of the Gaussian potential 
for (f)\ develops a vanishing moment. This is a first order phase transition |38| . 
Unlike generic first order phase transitions, it occurs at the point where the 
curvature of the potential first goes to zero. As a result, there is no energy 
barrier separating the two phases and no coexistence region of the ordered and 
disordered phases. 

These features of the phase transition are very sensitive to higher order 
corrections. In the following, we shall take into account the leading effect of 
higher order corrections. 

Since, as the temperature is raised from zero, the variable is the first 
to condense, we focus on its effective action. We have computed the relevant 
parts of the effective action up to three loop order. 

First of all, the ^-independent part of the effective action (what one obtains 
by putting cj) n = for vanishes up to order A 2 . 

The relevant parts which depend on (f) n have th.6 form 



Ai(r) = [l - 8r 3 - 3r 4 - 6r 2 ] - 24A [ln(r) r 2 (r 2 + l)(r + l) 4 ] - 

-3A 2 r 2 [ ln(r) 2 (68 r 10 + 352 r 9 + 904 r 8 + 1536 r 7 + 2256 r 6 + 
+3104 r 5 + 4120 r 4 + 2304 r 3 + 928 r 2 + 192 r + 16) - 
- ln(r) (27 r 10 + 152 r 9 + 390 r 8 + 640 r 7 + 915 r 6 + 1232 r 5 + 

+ 1748 r 4 + 1184 r 3 + 466 r 2 + 440 r + 102) ] + ...(24) 

A 2 (r) = i (1 + 8r 6 - 3r 8 - 6r 4 ) + ... (25) 





(26) 
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P 2 (r) = 3r 4 [ - ln(r) 2 ( 136r 12 + 512r n + 704r 10 - 1308r 8 - 1376r 7 + 1560r 6 
+6400 r 5 + 10896 r 4 + 8096 r 3 + 2136 r 2 + 1536 r + 240 ) + 
+ ln(r) ( 27 r 12 + 120 r 11 + 166 r 10 - 32 r 9 - 271 r 8 - 16 r 7 + 1044 r 6 + 
+2624 r 5 + 4036 r 4 + 3256 r 3 + 774 r 2 + 768 r + 944 ) ] + ... 

Eliminating 02 using its equation of motion, we obtain the effective action for 
01, in the large iV limit, and to order A 2 : 

j^S cS = Ai(r)|0i| 2 + A 2 (p 2 (r) - |0i| 4 + • • • (28) 

3.1 Phase transition 

As we raise the temperature from zero, the quadratic term in 0i in the effective 
action vanishes at the critical value of r, 



1 2° • I 
3 +A lF 



ln(3) - A 2 



2 2 -3 9 

which translates to the critical temperature 



11 -53 -3061 , , 13 2 -1867 1 , s5 
ln(3)+ 24 3§ ln(3)- 



+ 



/i 



121n(3) 



1 + A- 



2 b • 5 



3 1 



A 



23 • 19927 1765769 , , , 
+ -a ~ a ln(3) 



2 2 -3 7 



2 4 -3 8 



+ ... 



(29) 



(30) 



The zeroth order term in the critical temperature is the one found in 13 . 
The term of first order in A agrees with the result quoted in ref. |18j . 

Also, from (|28jl we see that the quartic term in <f>i is negative over the entire 
range < r < 1. This means that the phase transition is of first order. When 
r is just less than the critical rc = 1/3, the extremum of the effective action at 
01 = is only a local minimum. The effective action has a second zero when 



1 



Ai(r) 



i^i A 2 P 2 (r) -P 2 (r)/A 2 (r) 

Higher order terms in the effective action are individually small at this value 

of |0i I 2 when - , x tl^l,. , s « 1. We are further constrained by the fact 

Pi{r)— P( {r)/A2(r) J 

< A 2 << 1. The number 



P2(r)-Pj(r)/A 2 (r) 
that |0i | < 1. This requires that 



Ai(r) 
P 2 (r)-Pf(r)/A 2 (r) 



Ai(r) 
P 2 (r)-P?(r)/A 2 (r) 



is less than 0.10 in the range 0.2555 < r < 1/3 and is less 



than 0.001 in the range 0.3174 < r < 1/3. 
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If r is sufficiently close to r c , we can reliably say that the absolute minimum 
of the potential is not at <f>\ = but is elsewhere. This sets an upper bound on 
the transition temperature 

T CI it. < T c 



121n(3) 



The tunnelling barrier for bubble nucleation during the first order phase tran- 
sition is of order 1/ A 2 . 



4 Conclusions 

We have found that the phase transition in the weakly coupled plane wave 
matrix model is indeed of first order. As the temperature is raised from zero, the 
curvature contained in the quadratic term in the effective action still vanishes 
at some critical temperature. However, before that point is reached, when 
there is still an energy barrier between the two phases, the deconfined phase 
becomes the lower energy state. This is the generic behavior at a first order 
phase transition. In fact, this behavior is seen in other adjoint matrix models 
|41|-|43|. It is also the behavior that is seen in the collapse of Anti de Sitter 
space to a black hole, which is thought to be the analog of this phase transition 
in supergravity of a similar deconfinement in J\f = 4 supersymmetric Yang-Mills 
theory 

Our analysis does not allow us to compute the first order phase transition 
temperature accurately, only to deduce that it is of first order. It does, however, 
allow us to compute corrections to the Hagedorn temperature. This is the 
temperature at which, if the confining phase is superheated beyond where it 
is a global minimum of the free energy, it eventually becomes perturbatively 
unstable. It is just the place where the corrected curvature of the effective 
action vanishes. 

To conclude, we summarize our results: 

• We have computed the full free energy to two loop order: 
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, 2 -ioo P s _ 2W V ( (1-r 8 ) _ on (l-r 4 ) _ n (l-r 8 )(l-r 4 ) 
( r 8 + 4r 4 + 1)(r 4 _ 1)4 + [ cos/Mab + cos[3A bc + cos/3^ CQ ] 2r 4 (r 4 - l) 4 

1 s-iloi s~iloi 

U ab U bc °ca 

lg r 3 (r 4 - r 2 + l)(r 4 - l) 2 (r 2 + 1) [cos 0A ab + cos /M bc ] + r 6 (r 4 - l) 2 [2 + 2 cos (3 A, 



+32 



CabCbcCa 

r 3 (r 4 - l) 2 (r 2 + 1) [cos l3A ah + cos f3A bc ] + r 2 (r 8 - l)(r 4 - 1) cos f3A ca + (r 4 - l) 2 (r 8 + 1) 



C ab C bc C c 

Where we have: 



C2 = l-2r 4 cos/M afe + r 8 

C^ 2 = l-2r 2 cos/M afe + r 4 (31) 
C afe = l + 2r 3 cos/M^ + r 6 

This expression can be restated to exhibit its 4> n dependence: 



_L_ ^2-loops = 3A ^2 (j) n (j) m (j)_ rn _ n 



mn 

, (r 4 -r 2 + l)(r 2 + l) 2 3|ra|+4|m| _ Ifi ( r2 + 1 ) 2 / ^n+m 3|n|+3|m| 

+ i ° ( r 4 + r 2 + 1)(r 4 + 1 ) I l > 1D r 4 + r 2 + 1 ^ L ) 

! 30 ( r2 + !) 2 / 1 xn r 3|n|+2|m| _ 12 r 4|n|+2|m| _ 2Q r 2|n|+2|m| 

r 4 + r 2 + 1 

(- r 4 _ ]W 8 , 4 , / 2 , 10 _ 



r 4 + 1 



(r 4 + l) 2 



,4|n|+4|m| 
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(r 4 + l) 3 

(r 2 + l)(r 6 -r 4 + r 2 -l) 
(r 4 + r 2 + l) 2 



(r 4 + l)(r 4 + r 2 + l) 2 
(-l) m F mri (3,2) 

(32) 



Where we define the function F mn (a, b) in the following manner: 



F mn (a,b) = < 



m, n > or m, n < 

n < 0, m > —n or n > 0, m < —n (33) 

m < 0, m > —n or m > 0, m < — n 
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And where we have: 



(a,b) 



r ,a{2+n+m)+b 



„b(n+m)—an 



-b—2a+an 



1 



„2a+b 



+ 



„2a+b 



+ 



-2a—an+bn 



-2a—an+bn 



2a 



„— 2a— an+b(n+m) 



^—2a+an 
j ^*2ci 



+ 



(34) 



(a, 6) 



a(n+m)+6 



y,2a—an+b(n+m) ^, 
+ 



-b—an—bn 



+ 



1 _ r 2a+b ' ^ _ j,2a+fe 
j,— an+6(n+m) _|_ ^,—b—an ^ 

r b - 1 



-an 



—bn—b—an 

l-r b 



(35) 



—an+2a+bn ^,—b+an+2am ^,—an+bn 
_ j,2a+b ^ _ r 2a+b r b _ 

— an+£>(n+ni) ^— an+b(n+m) y(in+2om 
^-.o J ^.o _j_ ^-2a _j_ y-2a 

(36) 

In the zero temperature limit, r — > 0. We then see that the free energy is 
1 + 20 + 12 — 1 — 32 = 0, which is what is expected from supersymmetry. 

We obtain parts of the free energy to three loop order 

We check our three-loop computation by taking the zero temperature limit 
of it and finding that the free energy vanishes as the temperature is taken 
to zero. This is expected as a result of the supersymmetry of the model 
which is restored at zero temperature. 

We also find that the (^-independent part of the free energy vanishes to 
this order. Note that this is not the same as the zero temperature limit. 

We find the full shift in T c to order A 2 . It is given in eqn. (|3U|) . 

The form of the effective action that we find confirms the first order nature 
of the phase transition at weak coupling. 



Fl n {a : b)=r< n -^ +h 
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Appendix A: Three loops 

To do perturbation theory, we split the Euclidean action into three parts, a 
free action, an interaction action with three-point vertices and an interaction 
action with four-point vertices, 

S = / Tr -X 1 (-D 2 + ( M /6) 2 ) X 1 + -X~ a (-D 2 + (/i/3) 2 ) X~ a + 

J -(3/2 \ 2 2 

+^a ( _ £) + /i/4) ^ (37) 



rP/2 3 / i 
= / dTR2 Tr ifi-e^X" X b X~ c + ^ Ia af[X~ a , - 

J-r 



-\e a ^ aI g\j[X\^ J ] + l -e a ^ aI {g^) IJ [X\^ aJ ]\ (38) 



S 4 = - dr— Tr ([X\X j ] 2 + [X",X b ] 2 + 2[X°,X i } 2 ) (39) 
J 4 V / 



r /3/2 _ ^3 
-0/2 

respectively. 

The effective action Seff L4] is obtained by a perturbative expansion of the 
partition function. 



e -S eS [A] = J [dX][dip}e- s °- S3 - Si , S cS [A] = 1- < e- S:i ~ S4 > (40) 

where the bracket indicates the connected part of the expectation value in 
the free theory. Taylor expanding the exponential of interactions and retaining 
those terms which have non-zero contributions to two loop order gives 



n2 i if! 

U ab + 36 



det (-D ab ) det 8 {-D^ + f) 



+ 



+ [ < S 4 > ~ < si> 



1 

2 



\<Sl>-\< S 2 S A > +1 < Si > ) +...(41) 
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The first line on the right-hand-side is one- loop, the bracket after it contains 
two-loop and the next bracket the three-loop contributions. 

The appropriate two-loop Feynman diagrams are as follows, where a dotted 
line represents a fermion, and a solid line either flavour of scalar. 




At three loops there are many more diagrams. The following figures show the 
three loop diagrams and their associated zero temperature limits. The letter 
"P" denotes the circulation of the m = /x/3 scalar, while "Q" denotes the 
m = jtx/6 scalar. Where there is no indication, it is obvious from the allowed 
interactions dictated by the action. 

Cat's Eye Diagram 




Triple Bubble Diagram 
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Theta-Bubble Diagram 
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One can check that the sum of the above factors is zero, as garaunteed by 
SUSY. 

The vertices for these diagrams can be deduced from the interaction terms 
in the action. The propagators depend on temperature and are discussed in 
the appendix below. 

The technique that we use is to work in time, rather than momentum rep- 
resentation of the propagators. The integration over times is elementary and 
most of the work involves extracting the moments in r/ a = e l ^ Aa . Liberal use 
of computer algebra systems was used to unravel the expansion of diagrams in 
this object. 



Appendix B: Thermal Green functions 

The free field correlation function of the scalar field is 

[K b {T)Xl d {Q))=^5 ad 5 bc [ gi {r)] ab 
where the Green function is 
[<7i(r)L = (r 




(42) 



-T//3 



9(r)+ 



2 [rT T/ y(-r) 

where r\ = e tl3Aab = r] a /rj b and r = e _/3/1 / 12 . This Green function is defined 
on the interval t/(3 € (-1/2,1/2), with g\{l/2) = gi(-l/2) and it must be 
extended periodically to all real values r//3+integers. An alternative expression 
which is sometimes useful is 



9i{t) = — 



[r 2 ] |r|//3 + £ : 



41nr 



,2n 



V 



[r 2 ] T/f3 + ri n [r 



n=l 



The first term is the zero temperature Euclidean green function and the last 
two terms are the homogeneous solutions of the Euclidean wave equation which 
must be added to the first term in order to satisfy the periodic boundary con- 
ditions. Note that both green functions have a coefficient r/~ T /^. This factor 
cancels identically in all vacuum diagrams. It has been separated explicitly 
since some loop integrations are easier once it is canceled. 
Similarly, 
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ab 



[92 (r)} 



' ~D 2 ab + W3) 2 



|0) 



(43) 



81nr 



I r 41 r//3 , r/r 4 f 41 _r//3 



1 _ I r 4 [ 



+ 



1 — r/r 4 



0(t)+ 



l r 4 



+ 



1 - ±r 4 [ ; 
v 



+ 



1 — r\r L 



,41 -r//3 



0(-r) 



r 4lM/0 i 4n -n r 41 T //^ . n [ 



81nr 

The Fermion Green function is 
1 



„4i -r/jff 



71=1 



b/( T )L = ( 7 



-Z>a6 + (At/4) 



0) = -r? 



,31 r/i9 



1 + r 3 /ry 



)( r ) _ l r 3 0(-r) ) (44) 

7] 



It is antiperiodic 1/2) = —gf(l/2). 

Note that the expressions in l|42 |l .(|43j l are unchanged by the replacement 
r —s- 1/r. 
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